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Abstract 



Based on a high-resolution CO survey of Virgo spi- 
rals with the Nobeyama Millimeter-wave Array, we de- 
termined the dynamical centers using velocity fields, 
and derived position-velocity diagrams (PVDs) along 
the major axes of the galaxies across their dynamical 
centers. We applied a new iteration method to de- 
rive rotation curves (RCs), which reproduce the ob- 
served PVDs. The obtained high-accuracy RCs gen- 
erally show steep rise in the central ~ 100 to 200 pc 
regions, followed by flat rotation in the disk. We ap- 
plied a deconvolution method to calculate the surface- 
mass density (SMD) using the RCs based on two ex- 
treme assumptions that the mass distribution is either 
spherical or disk shaped. Both assumptions gave nearly 
identical results, agreeing with each other within a fac- 
tor of two at any radii. The SMD distributions re- 
vealed central massive cores with the peak SMDs of 
~ 10 4 — 10 5 Af ( g ) pc~ 2 and total masses of the order of 
~ 10 9 M ( g ) within a radius of 200 pc. Correlation analy- 
sis among the derived parameters shows that the cen- 
tral CO-line intensities are positively correlated with 
the central SMD, which suggests that the deeper is the 
gravitational potential, the higher is the molecular gas 
density in the nuclei regardless morphological types. 

Key words: galaxies: clusters of (Virgo) — galax- 
ies: rotation curve — galaxies: dynamics — Cluster: 
individual (Virgo) 



1 Introduction 



Rotation curves (RC) are one of the most basic tools 
to discuss the dynamics of galaxies and to derive the 



mass distribution. Extensive observations have been 
achieved in radio (HI: Bosma 1981a, b) and optical 
(Ha: Rubin et al. 1980, 1982, 1985, 1998; Mathewson 
et al. 1993, 1996). Outer flat RCs were noticed in the 
early studies (Rubin et al. 1980, 1982, 1985; Bosma 
1981a, b), and are often used to estimate the mass 
distributions in bulges, disks and halos (Rubin et al. 
1985; Kent 1986; Persic et al. 1996; Takamiya & Sofue 
2000). 

Detailed studies of central RCs of spiral galaxies have 
been started only recently (Sofue 1996, 1997; Sofue et 
al 1999; Bertola et al. 1998; Rubin et al 1997, 1999). In 
our studies of nearby galaxies, we have shown that the 
RCs generally rise steeply in the galaxy centers, and 
reach rotation velocities as high as ~ 100 — 300km s -1 
within radii of 100 - 200 pc. On the basis of these re- 
sults, we have conducted new observations of the cen- 
tral rotation properties of the Virgo CO rich spirals 
with high accuracy using the Nobeyama Millimeter- 
wave Array (NMA) in the 12 CO( J =1-0) line (Sofue 
et al. 2003a; Hereafter Paper I). 

In the decades, rotation curves have been derived 
from position-velocity diagrams (PVD) by optical (Ha 
, Nil) and radio (CO, HI) line observations (see, e.g., 
Sofue and Rubin 2001 for a review). There have been 
several methods to derive RCs from PVDs. The tradi- 
tional method is to trace intensity-weighted velocities 
(Warner et al. 1973), and the most popular method is 
to trace peak-flux ridges in PVDs (Rubin et al. 1985; 
Mathewson et al. 1992). These methods have given 
good results for outer regions of galaxies and for nearly 
face-on galaxies with sufficiently high spatial resolu- 
tions. 

Other methods are the terminal-velocity method 
used for our Galaxy (Clemens 1985), and the envelope- 
tracing method (Sofue 1996, 1997; Oiling 1996), which 
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traces the envelope velocities on PVDs and correct for 
the interstellar velocity dispersion and for instrumen- 
tal resolutions. This method gives better results for 
the central regions and for highly inclined and edge- 
on galaxies. We have applied this method to a number 
of intermediately inclined galaxies to obtain central-to- 
outer RCs (Sofue 1996; Sofue et al. 1999). 

These various methods have been the most practi- 
cal ways to derive RCs from PVDs. However, they 
do not necessarily give good results for the central re- 
gions, where the PVDs show large intrinsic velocity 
width and high interstellar velocity dispersion. The 
central PVDs also suffer strongly from the smearing 
effect due to finite beam sizes and finite slit widths. 
In fact, it is known that RCs derived by these meth- 
ods cannot reproduce the observed PVDs in the in- 
nermost regions of galaxies, when they are convolved 
with the observed intensity distributions. In order to 
resolve these problems, we have developed a new it- 
eration method (Takamiya and Sofue 2002), which we 
describe in detail and apply to our new data in Section 
2. 

Not only the methods, but also the properties of the 
used lines and observational conditions have limited 
the accuracy of the central RCs. Resolutions in HI 
observations were usually tens of arcseconds and are 
insufficient to resolve the central disks, and also the HI 
gas is often deficient in the centers. Optical lines are 
severely absorbed by the central dusty disks, particu- 
larly in cases of edge-on galaxies. Optical spectra are 
severely contaminated by the bright continuum emis- 
sion of the bulge and nucleus. The Ha emission is 
often superposed by a broad absorption line due to the 
Balmer wing, and is often superposed by an extremely 
broad emission line due to Seyfert activity. 

In this context, we have stressed that the CO line is 
an ideal tool to investigate the central kinematics and 
RCs because of the following advantages (Sofue 1996, 
1997; Sofue et al. 1999): (1) CO gas is strongly con- 
centrated in the central region. (2) The CO gas has 
much smaller intrinsic velocity dispersion than stars, 
so that the velocities manifest the rotation. (3) Even 
though each molecular cloud is optically thick for the 
CO line, the high-velocity rotation and velocity dis- 
persion among clouds make the line practically trans- 
parent in so far as the galactic-scale dynamics is con- 
cerned. (4) High-angular resolution can be achieved 
by mm-wave interferometry. (5) High-velocity resolu- 
tion can be achieved by radio spectrometers combined 



with heterodyne receivers. (6) High sensitivity can be 
achieved by large- aperture telescopes and interferome- 
ters. 

The Virgo high-resolution CO survey (ViCS) has 
been conducted in order to obtain a high angular- and 
spectral-resolution database for a large number of CO- 
bright Virgo Cluster spirals in the 12 CO (J = 1 — 0) line 
using the Nobeyama Millimeter-wave Array (NMA)of 
NRO 1 . The major scientific motivation was to inves- 
tigate detailed central kinematics of the galaxies, par- 
ticularly the innermost RCs from analyses of PVDs 
across the nuclei, and to derive the central mass distri- 
butions. The advantage to observe the Virgo member 
galaxies is the identical distance of 16.1 Mpc based on 
the Cepheid calibration (Ferrarese et al. 1996). So, 1 
arcsec corresponds to 78 pc, and 1 kpc corresponds to 
12".8. 

The observations were obtained in the course of a 
long-term project during the winter seasons from 1999 
December through 2002 April in the AB (longest base 
lines), C (medium) and D (short base lines) array con- 
figurations of the NMA. The resultant spatial resolu- 
tions were 2" to 4", and velocity resolutions were 10 
to 20 km s _1 . The detailed description of the survey 
together with the observational parameters and reduc- 
tion procedure are given in Paper I (Sofue et al. 2003). 

In this paper we use the observed PVDs to derive 
RCs by applying the iteration method, and calculate 
radial distributions of the surface-mass density (SMD) 
by applying a SMD-deconvolution method described in 
Takamiya and Sofue (2000). 

2 Rotation Curves using the It- 
eration Method 

2.1 Iteration Method 

RCs have been derived from observed PVDs, and thus, 
should be able to reproduce the original PVDs when 
convolved with intensity profiles. The iteration method 
(Takamiya & Sofue 2002) determines a RC so that it 
becomes consistent with observations on PVD-basis. 
We describe the detailed procedure below. The algo- 

1 NRO (Nobeyama Radio Observatory) is a branch of the Na- 
tional Astronomical Observatory, operated by the Ministry of 
Education, Culture, Sports, Science and Technology, Japan. 
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rithm is illustrated in a flow chart in figure 1. We 
assume axisymmetry of a galactic disk throughout the 
iteration. For a simple explanation, we describe below 
an observed PVD as PV b s , a model PVD and RC at 
the i-th iteration as PV \ and Vi, respectively. 

- Fig. 1 - 

Before the iteration, we define a reference RC (V rc f) 
and a three-dimensional intensity distribution (/ re f ) in 
a galaxy disk. V rc f is defined as a radial velocity profile 
traced at the 20%-level envelope of the peak flux in the 
PVobs- The intensity, / ro f, is assumed to have a form: 

I le f(r,z) = AS(r)scch 2 (z//i z ), (1) 

where S(r) is obtained by integrating PV^bs along 
the velocity at each radius, sech 2 (z/h z ) is a model z- 
direction structure, and A is a normalization constant. 
We assume the scale hight of h z = fOO pc, and the 
intrinsic interstellar velocity dispersion is isotropic and 
is a v — 10 km s _1 in the gas disk. An initial guess 
of RC (Vb) is set to be Vq = Kef- Then the iteration 
proceeds as follows. 

Step 1: Constructing a model galaxy by convolving 
Vi and iref with taking a v into account. 

Step 2: Virtually observing the model galaxy to ob- 
tain PVi in the same conditions (slit-width, resolution, 
inclination, etc.) that PV bs was observed. 

Step 3: Determining a temporal RC (K) by tracing 
the 20%-level envelope of the peak flux in PV^. 

Step 4: Calculating the difference 

AVi = V ie{ -V t , (2) 

which should be close to the difference between K and 
the real RC that we want to obtain. 

Step 5: Applying an iteration termination criteria, 

|AK|<<7 tol , (3) 

where er to i is an arbitrary tolerance velocity; we set it 
as cr to i = 20km s _1 . If the criterion is not satisfied, wc 
go to Step 6. If the criterion is satisfied, the iteration 
terminates, and the resultant RC is Vi. 

Step 6: Calculating the K+i for the next iteration 
step (i + 1) as 

K+i = Vi + AVi, (4) 



and we go back to Step 1 . 

We made modification to the original method of 
Takamiya & Sofue (2002) in treating equation (2). In 
the original method, only the radial velocity profiles at 
the 20%-level envelope in PV b s and PVi were consid- 
ered. Wc further take the 20, 40, 60, 80, 100%-level 
envelopes (K r ' e/ and K/; I = 20, 40, 60, 80, 100%), and 
replace the equation (2) with 

AV = ^>-(Ki f -Vi), (5) 
i 

where the weights gi satisfy J2i9i — 1- We usually 
set <?i = .92 = ••• = 9u while some g\ s are occasionally 
set to be zero in cases that those V 1 are ill-determined 
owing to noises. This modification avoids instability 
occurred in the iteration that comes from fluctuations 
on the 20%-level envelopes. When g t = 1 for I = 20% 
and for the others, equations (2) and (5) become 
identical. 



2.2 Center Positions 

The center positions are adopted from Paper I, where 
we applied the task GAL in the AIPS reduction pack- 
age, which fits the velocity field (the first moment map 
of the RA-Dec- Velocity cube), by a Brandt rotation 
curve (Brandt and Scheer 1965). Table 1 lists the thus 
obtained dynamical centers and systemic velocities. In 
most cases the observed centers of the fields of view co- 
incided well with the observed centers within one arc- 
second. However, some deviation was found for NGC 
4254 and NGC 4402, and we reset the centers of these 
galaxies to the determined dynamical centers. 

The GAL fit gave also the inclination and position 
angle of the inner disk. However, these quantities are 
dependent on the streaming motions superposed on the 
circular rotation, which are not negligible in the cen- 
tral regions. Hence, we adopted the inclination and 
position angles of the whole galaxies determined from 
optical images of wider disks. In table 1 we list the 
adopted inclination angles taken from Paper I and the 
literature therein. We adopt these values for construct- 
ing PVDs and rotation curves. 
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Table 1: Adopted center positions, systemic velocities, and parameters for RC-PVD iteration 



Galaxy RA (J2000) Dec. (J2000) V sys lsr Inclination Beam size Slit width for PVD 
(km s- 1 ) (°) (") n 



NGC 


4192 


12hl3m48. 


29s 


+14d54m01. 


,9s 


-127.3 


74.0 


2, 


,0 


3. 


.0 


NGC 


4254 


12hl8m49. 


,61s 


+ 14d24m59. 


6s 


2410.8 


28.0/42.0 


2, 


.5 


3, 


.0 


NGC 


4303 


12hl9m21. 


67s 


+04d45m03. 


,7s 


1557.9 


25.0 


2, 


.5 


3, 


.0 


NGC 


4402 


12h26m07. 


,44s 


+ 13d06m44. 


,7s 


216.4 


75.0 


2. 


.5 


3, 


.0 


NGC 


4419 


12h26m56. 


,40s 


+ 15d02m50. 


,2s 


-188.6 


67.0 


3, 


,0 


3 


.0 


NGC 


4501 


12h31m59. 


,12s 


+ 14d25ml3. 


,3s 


2261.1 


58.0 


4, 


.5 


3, 


,0 


NGC 


4535 


12h34m20. 


,35s 


+08dllm52. 


,2s 


1966.3 


43.0 


3 


,0 


3, 


.0 


NGC 


4536 


12h34m27. 


,08s 


+02dllml7. 


,0s 


1801.7 


67.0 


2, 


.2 


3 


,0 


NGC 


4548 


12h37m57. 


,47s 


+ 14dl3ml8. 


,0s 


464.9 


37.0 


2, 


,3 


3, 


.0 


NGC 


4569 


12h36m49. 


,79s 


+ 13d09m48. 


,7s 


-195.8 


63.0 


4, 


.0 


3, 


,0 


NGC 


4654 


12h41m25. 


,59s 


+ 13d24m03. 


,2s 


1051.4 


52.0 


5, 


,0 


3 


.0 


NGC 


4689 


12h50ml5. 


,86s 


+ 13d29m27. 


,4s 


1614.8 


30.0 


5, 


,0 


3 


.0 



2.3 Position- Velocity Diagrams 

In order to obtain PVDs of the observed galaxies, we 
first rotated the data cube in (RA, Dec, Velocity) to 
(X, Y, V) , so that the first axis X represents the dis- 
tance in the major axis direction, the second axis Y the 
distance along the minor axis, and the third axis V the 
radial velocity. We then transposed it to a cube with 
the axes (X, V, Y). Here, each (X, V) plane represents 
a PVD sliced at a distance Y from the minor axis. 

The PVD along the major axis was obtained as that 
at Y = 0. In order to increase the signal-to-noise ratio, 
we averaged several diagrams within Y = ±1".5 from 
the major axis, or those within a slit of 3" wide across 
the dynamical center. However, the beam sizes were 
comparable to or greater than the slit width as shown 
in table 1. The effective spatial resolutions in the PVDs 
are, therefore, (beam width 2 + slit width 2 ) 1 / 2 , which 
are typically a few seconds of arc (1" = 78 pc, and 1 
kpc = 12". 8). The upper diagrams of figure 2 show the 
thus constructed observed PVDs along the major axes 
across the dynamical centers after correcting for the 
inclinations of galaxy disks. The intensities are grey- 
scaled, so that the maximum intensity is in black, and 
the minimum in white. 

- Fig. 2 - 



2.4 Rotation Curves 



We applied the iteration method to the observed PVDs, 
and obtained rotation curves for 12 galaxies listed in 
table 1. Both sides of each PVD along the major 
axis were fitted independently. The obtained RCs are 
shown in the lower diagrams of figure 2, where we 
show the reproduced PVDs calculated by convolving 
the RCs with the observed intensity distributions. In 
most cases, the convolved PVDs well mimic the obser- 
vations. We then averaged the absolute values of the 
RCs from both sides for each galaxy, and obtained ro- 
tation curves as a function of the radius. We combined 
the thus obtained CO-line RCs with outer optical RCs 
from Rubin et al. (1999). We also used outer RCs from 
Sofue et al. (1999) for NGC 4303 and NGC 4569 for 
larger radii, where both optical and HI data have been 
combined. When there are two or more observations 
from the literature, we took the average. Since the 
inclination angles used in the literature are not neces- 
sarily the same as those used in our analysis, we have 
corrected them to the same inclination angles as in ta- 
ble 1. Figure 3a shows the combined RCs of the Virgo 
galaxies, where the inner CO RCs from the present 
observations are drawn by thick lines. In figure 3b we 
plot all RCs in one diagram in order to see their general 
property. 
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The resultant RCs in figure 3 indicate very steep rise 
in the central few arcsec regions, or in radii 100 to 200 
pc. The steep rise near the center is remarkable when 
we compare with those from Ha and HI observations. 
Note, however, that fluctuations of about ±10 km s _1 
superposed on the RCs are not real, which arise mainly 
from patchy distributions of emission in the PVDs. 

- Fig. 3 - 



3 Mass Distributions 

Radial distribution of the surface-mass density (SMD), 
a(R), can be calculated directly from a rotation curve 
(Takamiya and Sofue 2000, and the literature therein). 
We may assume that the 'true' mass distribution in a 
real disk galaxy will be in between two extreme cases: 
spherical and axisymmetric flat-disk distributions, as 
discussed in detail in Takamiya and Sofue (2000). 

3.1 Spherical Mass Distribution 

If we assume that the mass distribution is spherically 
symmetric, the SMD <rs(R) at a radius R is calculated 

by, 

CO 

a s (R) = 2 J p{r)dz, (6) 
o 

oo 

= / , 1 dM } r) d, (7) 
2tt J rVr 2 - R 2 dr 

R 

Here, R, r and z are related by r = \/R 2 + z 2 , p(r) the 
spatial density, and M(r) is the mass inside radius r 
given by 

M(r) = (8) 

where V(r) is the rotation velocity at r given by the 
rotation curve. 

Note that the equation gives underestimated mass at 
R <~ -R m ax because of an edge-effect that the integra- 
tion is done only up to a maximum radius of observa- 
tion, i? max , instead of infinitely large R. However, this 
edge effect is negligible for the calculated result within 
a radius slightly smaller than i? max , because the factor 
(•s/r 2 — R 2 ) rapidly diminishes as r gets away from 
R. 



3.2 Flat-Disk Mass Distribution 



If we assume that the galaxy comprises an infinitely 
thin flat disk, the SMD, ar>{R), can be calculated by 
solving the Poisson's equation: 




where K is the complete elliptic integral and becomes 
very large when x ~ R (Binney & Tremaine 1987). 



When we calculate ctd (R) , we must take into account 
the following points. First, the equation is subject to 
the boundary condition, V(0) = V(oo) = 0. Because 
we have only a few data points for the very central re- 
gion, where the assumption of V(0) = may not apply 
from the data, the calculated result near the nucleus 
would not be reliable. Hence, in so far as the central 
<~ 100 pc is concerned, the spherical assumption will 
give more reliable SMD. Second, the upper boundary 
of the integration of the second term is i? ma x, where 
V(R) is finite, instead of R =infinity, where the equa- 
tion assumed zero velocity. This leads to overestima- 
tion of <jt>(R) at R ~ -Rmax, and the SMD near the 
outer edges will give an upper limit. 

3.3 Distributions of the Surface-Mass 
Density 

Figure 4a shows the thus calculated SMDs for indi- 
vidual galaxies, where thick dotted lines represent the 
results for flat-disk assumption, and thin full lines 
for spherical assumption. Since the central masses of 
galaxies are considered to be dominated by spheroidal 
components, we may rely more on the SMDs from 
spherical assumption. On the other hand SMDs from 
flat-disk assumption give better results for outer re- 
gions, where the disk component dominates. So, in 
the following, we discuss the SMD's properties in the 
central regions inside r ~ 1 kpc based on the SMDs 
from spherical assumption, and those for outer disks 
at r >^ 1 kpc on those from flat-disk assumption. 
In figure figure 4b we show all the SMD profiles from 
flat-disk assumption in one diagram in order to see 
their general property. Figure 4c shows all SMDs from 
spherical assumption for the central 2 kpc. 
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- Fig. 4 - 

In spite of the variety of galaxy type and morphology 
of molecular gas distributions (Sofuc et al. 2003), the 
SMDs appear to show a basic, principal structure as 
the following. 

(1) Central massive core: The mass distribution is 
highly peaked at the center, showing a core component 
with scale radius of about 100 to 200 pc. Since the es- 
timated radius is comparable to or smaller than the 
beam size, this component is likely not resolved, but 
may have much smaller radius. The calculated central 
SMDs amount to ~ 10 4 — 10 5 M,gpc~ 2 , which is, how- 
ever, lower limit, because the core components are not 
sufficiently resolved. The dynamical mass of this com- 
ponent within 100 to 200 pc radius is of the order of 
~ 10 9 Af®. In table 2 we list the determined SMDs at 
the centers and dynamical masses within radii 200 and 
500 pc. 

(2) Bulge component: The SMD is then followed 
by a more gradually decaying profile at r = 0.2 — 2 
kpc with scale radius of 500 pc. This component is 
due to the central bulge. A typical case of distinct 
components is seen for NGC 4536. However, in some 
galaxies the bulge and core component are not clearly 
distinguishable because of the insufficient resolution. 

(3) Disk component: The radial profiles at radii 
greater at r ~ 2 to ~ 7 kpc are approximately ex- 
pressed by an exponentially decaying function of scale 
radius of a few kpc. 

(4) Massive halo: At radii r >~ 7- 8 kpc, the SMD 
profile becomes more gently decreasing, indicating an 
excess over the exponential disk due to the massive 
halo. This is clearly seen in galaxies with RCs covering 
larger radii. 



3.4 Correlation among dynamical 
properties 

In table 2 we list the morphological types of the an- 
alyzed galaxies, their total B magnitude taken from 
NED, and the peak SMD at the center, S c . We also list 
dynamical masses within 200 pc and 500 pc radii from 
the dynamical centers, M 2 oo and M 50 o, which were cal- 
culated by 

M R = RV? ot /G, (10) 



where V ro t is the rotation velocity given by the RCs, 
and G the gravitational constant. 

SMD vs B magnitude: Figure 5a shows a plot of the 
center values of SMD (column 4 of table 2) against total 
B magnitude (column 3) . There is a trend that brighter 
galaxies have higher SMD. This is consistent with the 
fact that brighter galaxies have brighter bulges, as 
shown below for the bulge mass and central SMD. 

Mass vs B magnitude Figure 5b is a plot of dynamical 
masses within radii 200 pc (column 5 of table 2; filled 
circles) and 500 pc (column 6; diamonds) against total 
B magnitude. The trend is the same as for SMD vs B. 

SMD vs Dynamical mass: Figure 5c plots the cen- 
ter SMD (column 4 of table 2) against the dynamical 
masses within radii 200 (column 5; filled circles) and 
500 pc (column 6; diamonds). There is a tight corre- 
lation between the mass and SMD, which are related 
as 

S c cx M 2 ° 8 o, (11) 

and 

S c cx M 500 . (12) 

This correlation is naturally understood, if the core 
radii of the central mass concentrations are not largely 
different among the galaxies. 

Since the rotation curves are nearly flat at radii 0.2 
to 1 kpc in most galaxies, these dynamical masses and 
center SMD are positively correlated with the bulge 
mass, if the bulge radius is of the order of several hun- 
dred pc. Figure 5d plots the dynamical mass within 
200 pc radius against that within 500 pc from the same 
table, which shows a linear correlation. The correla- 
tion between the center SMD and central mass within 
200 pc with the surrounding dynamical mass would be 
somehow related to the correlation found for the masses 
of bulges and the central black holes (Kormendy and 
Richstone 1992). 

— Fig. 5 a,b,c,d — 

3.5 Correlation between CO concentra- 
tion and dynamical properties 

The seventh column of table 2 lists the integrated CO- 
line intensities at the centers determined by ellipse fit- 
ting, Ico.c using the observed first moment maps, as 
taken from Paper I. The last column lists the peak CO 
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Table 2: Derived mass parameters and CO intensities. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Galaxy 


Morphology 


Btot 


Center SMD 


M 200 


M 500 


Ico.c cosi 


Ico, P cosi 






[mag] 


[lO^pc- 2 ] 




[10 9 M®] 


[K km s- 1 ] 


[K km s- 1 ] 


NGC 4192 


SAB(s)ab 


10.95 


3.4 


1.16 


6.42 


73 


96 


NGC 4254 


SA(s)c 


10.44 


1.4 


0.32 


1.64 


60 


82 


NGC 4303 


SAB(rs)bc 


10.18 


1.9 


0.69 


1.23 


160 


395 


NGC 4402 


Sb 


12.55 


0.5 


0.14 


1.03 


5 


58 


NGC 4419 


SB(s)a 


12.08 


1.3 


0.34 


1.67 


100 


229 


NGC 4501 


SA(rs)b 


10.36 


2.0 


1.06 


4.51 


88 


111 


NGC 4535 


SAB(s)c 


10.59 


2.0 


0.57 


3.36 


200 


275 


NGC 4536 


SAB(rs)bc 


11.16 


3.0 


1.10 


3.72 


171 


224 


NGC 4548 


SBb(rs) 


10.96 


4.2 


1.92 


6.92 


38 


98 


NGC 4569 


SAB(rs)ab 


10.26 


0.9 


0.15 


2.31 


162 


266 


NGC 4654 


SAB(rs)cd 


11.10 


0.35 


0.06 


2.06 


10 


27 


NGC 4689 


SA(rs)bc 


11.60 


0.75 


2.98 


0.13 


12 


21 



(1) NGC number; (2) Morphological type from RC3; (3) Total B magnitude from NASA Extragalactic Database 
(NED); (4) Surface mass density (SMD) at the center; (5) (6) Dynamical mass calculated by spherical assump- 
tion. M200 and M500 are masses within 200 and 500 pc radii, respectively. (7) CO intensity at the center 
calculated by ellipse fitting, taken from Paper I corrected for inclination. (8) Peak CO intensities in the ob- 
served first-moment maps, taken from Paper I corrected for inclination. Note that the CO peaks are near to 
the centers, but do not necessarily coincide with the center positions. 
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intensities in the maps, Ico,p- The CO intensities arc 
corrected for the inclination. Note that the CO inten- 
sity peaks do not necessarily coincide with the dynam- 
ical centers, which are, however, located near to the 
centers. 

CO intensity vs B magnitude: Figure 6a plots the 
central CO intensity (column 7 of table 2; circles) and 
peak CO intensity (column 8; triangles) against total 
B magnitude (column 3). There is a weak correlation 
between B luminosity and central and peak CO inten- 
sities: brighter galaxies have higher CO intensities. 

CO intensity vs SMD: An interesting correlation was 
found in figure 6b, where the center CO intensity (col- 
umn 7 of table 2; circles) and peak CO intensity (col- 
umn 8; triangles) are plotted against the center SMD 
(column 4). There is a clear correlation that the higher 
is the SMD, the higher is the peak and central CO in- 
tensities. This fact indicates that the deeper is the 
central gravitational potential, the more strongly the 
CO gas is concentrated in the nuclei. 

- Fig. 6 a,b - 

4 Description of Individual 
Galaxies 

We describe properties of individual galaxies seen in 
their PVDs, rotation curves, and mass distributions. 
More general description of individual galaxies, includ- 
ing CO intensity distributions and velocity fields, are 
presented in Paper I. 



4.2 NGC 4254 

The velocity field shows a regular spider pattern, indi- 
cating a circular rotation of the disk. The PVD and 
RC show a sharp rise of rotation velocity in the central 
few arcseconds, reaching a shoulder-like step at 100 km 
s _1 , and then the velocity increases gradually to 150 to 
200 km s _1 , depending on the inclination angle. Note 
that there are two possibilities of inclination angles, 
cither 28°infcrrcd from molecular gas disk or 42°from 
K-band images (Sofue et al. 2003b). The SMD shows 
clearly a massive core with a sharp peak and a bulge 
component. They are surrounded by a disk compo- 
nent, while the disk part is uncertain for the limited 
coverage in the present data. 



4.3 NGC 4303 

The intensity distribution is classified as twin-peaks, 
and shows two prominent off-set bar ridges (Koda et 
al. 2003). Due to the non-circular motion, the PVD 
has a wider velocity dispersion compared to the other 
galaxies. The RC indicates steep rise of rotation within 
a few arcseconds, but it is wavy at r < 3 kpc radius. 
Such wavy behavior of RC is typical for a barred galaxy 
(Sofue et al. 1999). Accordingly, the SMD behaves 
peculiar at 1 to 4 kpc, and no firm conclusion about the 
inner mass distribution can be obtained. However, the 
velocity field in the very central region within radius 
300 pc shows a regular spider pattern, where the steep 
rise of rotation is observed. Hence, the central massive 
component shown in the SMD would be real, as well 
as that at r > 5kpc. 



4.1 NGC 4192 

A bright CO intensity peak is observed at the center, 
and the galaxy is classified as a "central/single-peak" 
type (Paper I). The rotation characteristics appears 
normal with high velocity rotation at ~ 500 pc. The 
rotation velocity, then, declines to a minimum at 3-4 
kpc according to the outer rotation curve taken from 
the literature (Sofue et al. 1999). It is then followed 
by a flat rotation, which continues to the outer edge. 
The SMD shows that the central core and bulge com- 
pose a single component, which is surrounded by an 
exponential disk. 



4.4 NGC 4402 

The CO intensity distribution shows a high density nu- 
clear disk of radius ~ 10", which is surrounded by 
a more broadly distributed molecular disk with arm 
structures. The PVD and RC shows steeply rising nu- 
clear rotation, followed by a gradually rising outer disk 
rotation. The SMD indicates a sharp peak at the cen- 
ter, suggesting a massive core with small scale radius, 
while the present resolution is not sufficient. It is then 
surrounded by a bulge and disk components. 
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4.5 NGC 4419 



4.10 NGC 4569 



The CO gas is strongly concentrated in the central 5" 
radius disk, showing a single-peak surrounded by outer 
disk. The central disk has steeply rising rotation, fol- 
lowed by a gradual rise of the disk component. The 
SMD comprises a central core plus bulge component, 
while the disk part does not show up because of the 
insufficient radius coverage. 

4.6 NGC 4501 

A compact central peak and spiral arms are promi- 
nent in the intensity map. The velocity field is normal, 
showing a typical spider pattern. The PVD shows a 
central component superposed by a disk/arm compo- 
nent. The RC indicates a sharp central rise followed 
by a gradually rising part and flat outer rotation. The 
SMD indicates a central core plus bulge and disk. 

4.7 NGC 4535 

A high CO gas concentration is found at the center, 
and the velocity field is normal, showing a spider pat- 
tern. The PVD comprises a single ridge, and the RC 
rises steeply near the center, reaching a shoulder like 
maximum followed by a wavy frat part in the disk. 
SMD indicates a core, bulge and disk. 

4.8 NGC 4536 

This is also a single- CO peak galaxy with a typical 
spider velocity field. The PVD comprises two ridges, 
a central steeply rising component and a rigid-body 
like disk component. The RC rises extremely steeply 
in the central 100 pc, reaching a maximum followed 
by a flat RC. SMD shows clearly distinguishable three 
components of the core, bulge and disk. 

4.9 NGC 4548 

CO gas is concentrated in the central few arcsec region, 
and the intensity is weak. However, the PVD and RC 
show very high rotation in the central several hundred 
parsecs. Accordingly, the SMD shows a high-density 
central core, bulge and disk components. 



CO gas is distributed in an elliptical ring of radius 500 
pc (Nakanishi ct al. 2003), and the velocity field is com- 
plex, having anomalously high non-circular rotation. 
The PVD is complicated with two bright ridges having 
forbidden velocities. Our iteration program averaged 
the observed PVD, and obtained an RC, which gives 
only averaged rotation velocities. The SMD shows a 
bulge component and disk, while its detail in the cen- 
tral 1 kpc is not reliable because of the uncertainty in 
the RC determination. 



4.11 NGC 4654 

CO distribution is highly distorted and lopsided. The 
PVD shows rigid-body ridge, and the RC is also grad- 
ually rising at r < 2kpc. Accordingly, the SMD com- 
prises only disk and bulge components with slight en- 
hancement due to a possible core component near the 
center. 



4.12 NGC 4689 

CO gas is broadly distributed with amorphous mor- 
phology, but the velocity field is normal, showing a 
spider diagram. The CO intensity is weak and its dis- 
tribution is patchy. PVD has a central ridge and ex- 
tended flat part. RC shows a sudden rise in the central 
~ 100 pc to a velocity of ~ 70 km s _1 , and then rises 
gradually to a flat part. No clear compact component 
is seen in the SMD, while bulge and disk components 
show up. 



5 Discussion 

We have analyzed the high-resolution CO-line survey of 
Virgo spirals obtained with the Nobcyama Millimeter- 
wave Array (Paper 1) in order to derive dynamical pa- 
rameters of the galaxies. The study of Virgo galaxies 
provided us with a unique opportunity of using a very 
accurate, unambiguous Cepheid distance of 16.1 Mpc 
from the HST observations (Ferrarese et al. 1996). 

We determined the dynamical centers using the ob- 
served CO-line velocity fields obtained by the high res- 
olution CO survey of CO-rich Virgo spirals. The high 
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concentration of molecular gas in the centers and high- 
spatial and velocity resolutions enabled us to determine 
the dynamical centers very accurately. From the data 
cube, we constructed position-velocity diagrams along 
the major axes across the dynamical centers. 

By applying the new iteration method to derive 
exact rotation curves from PVDs as proposed by 
Takamiya and Sofuc (2002), we determined rotation 
curves. Thus obtained RCs can well reproduce the 
observed PVDs by convolving with the observed CO 
intensity distributions. We combined these RCs with 
outer disk RCs from the literature. Except for a few 
cases, the RCs rises steeply near the center within the 
spatial resolution. Attaining a central peak or a shoul- 
der after the rapid rise, the RCs are followed by the 
gradually varying disk RCs which are in general flat 
till the outer edge. Some galaxies like NGC 4303 and 
NGC 4659 with strong non-circular velocity field due 
to a bar showed wavy behavior in the central ~ 1 kpc. 

We, then, calculated the distributions of surface- 
mass density SMD as a function of radius on the two 
extreme assumptions of spherical symmetry and flat- 
disk using the method developed by Takamiya and So- 
fue (2000). Both the spherical and disk results agreed 
within a factor of two, while the spherical assumption 
give more reliable results for the innermost regions, 
where bulge and core component are dominant. On 
the other hand, the flat-disk assumption give better 
results for the disk and outer regions where the disk 
component dominates. In the present study the halo 
component was not thoroughly detected, because we 
have concentrated in the analysis of the innermost RCs. 

The SMDs generally have strong central condensa- 
tion of mass of scale radii of 100 to 200 pc, which we 
call the massive central core. The core radii are signif- 
icantly smaller than those of bulges, whose scale radii 
are usually several hundred parsecs to one kpc. Dy- 
namical masses of the massive cores within a radius 200 
pc are of the order of 10 9 M^. Since the present spatial 
resolutions (2 — 4", 160 -300 pc) are comparable to the 
supposed scale radii, we may have not fully resolved 
the core components. Therefore, the scale radii may 
still be smaller, and hence, the center SMD derived by 
the present analysis would give lower limits. However, 
the fact that the center SMD is roughly proportional 
to the dynamical mass within 200 pc radius (figure 5) 
suggests that the scale radii of the cores would not be 
largely scattered among the galaxies. 



The present results about the SMD distribu- 
tions confirm our earlier results about massive cores 
(Takamiya and Sofuc 2000; Sofue ct al. 2001; Koda et 
al. 2002). Koda ct al. (2002) argued in detail that the 
high- velocity rotation observed in the central region of 
NGC 3079 indeed manifests the central massive core 
with deep gravitational potential. 

An interesting positive correlation was found be- 
tween the CO gas concentration in the center and the 
central SMD. The observed galaxies are not dominated 
by barred galaxies, and moreover, the CO concentra- 
tion is not dependent on galaxy types. This fact implies 
that the often mentioned inflow mechanism of molec- 
ular gas due to bar-driven angular momentum loss is 
not absolutely necessary. In fact, our analysis in Paper 
I showed that the highest CO densities are found more 
in centrally single CO peaked galaxies than in twin- 
peaked galaxies indicative of bar driven-inflow. The 
Ico - SMD correlation indicates simply that the deeper 
is the central potential, the stronger is the gas concen- 
tration. 

This is consistent with our recent results about the 
extremely high density molecular core in NGC 3079, 
which was argued to be the consequence of a very deep 
potential, where the dense gas is maintained gravita- 
tionally stable due to high velocity shears, and is sup- 
pressed from being fragmented to form stars (Sofuc 
et al. 2001; Koda et al. 2002). The mechanism of 
molecular gas concentration and its maintenance with- 
out suffering from being exhausted by star formation in 
a deep gravitational potential would be an interesting 
subject for the future. 
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Figure Captions 

PS figures are avaiable from http://www.ioa.s.u- 
t okyo . ac . j p /radio / virgo / 

Fig. 1. Algorithm of the iteration method to derive an 
exact rotation curve from an observed position- velocity 
diagram (Takamiya and Sofue 2002). 

Fig. 2. Upper panels: Observed PVDs along the ma- 
jor axes after correcting for the inclination of galaxy 
disks. Lower panels: Obtained rotation curves by the 
iteration method, and reproduced PVDs by convolving 
the RCs with observed intensity distributions. Here, 
1" = 78 pc and 1 kpc =12". 8 for the adopted Virgo 
distance of 16.1 Mpc. 

Fig. 3: (a) Rotation curves of the Virgo spiral galaxies 
obtained by the iteration method, combined with outer 
RCs from the literature (Rubin ct al. 1999; Sofue et 
al. 1999), for which inclination angles are corrected to 
the same values as in Table 1. The CO RCs from the 
present observations are drawn by thick lines. 

(b) All rotation curves in one diagram. 

Fig. 4. (a) Radial profiles of the SMD (surface-mass 
density) of the Virgo spirals. Full lines are the results 
from spherical symmetry assumption, which are more 
reliable in inner region. Dashed lines are the result 
from flat-disk assumption, which are more reliable for 
the disk and outer region. 

(b) All SMD profiles from flat-disk assumption plotted 
in one diagram. 

(c) All SMD proficls from spherical assumption for the 
central 2kpc region plotted in one diagram. 

Fig. 5 (a) SMD at the centers plotted against total B 
magnitude. 

(b) Dynamical masses within r = 200 pc (dots) and 
500 pc (diamonds) against total B magnitude. 

(c) SMD against dynamical masses within r = 200 
(dots) and 500 pc (diamonds). 
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(d) Dynamical mass within r = 200 pc against mass 
within r = 500 pc. 

Fig. 6. (a) Peak (triangles) and central (circles) CO in- 
tensities corrected for inclination angle plotted against 
total B magnitude. 

(b) Peak (triangles) and central (circles) CO intensities 
against the central SMD. 
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